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I. INTRODUCTION 



Heavy quarkonium plays an important role in establishing and understanding quantum 
chromodynamics (QCD), the fundamental theory of strong interactions. Due to the existence 
of several energy scales involved with these systems, heavy quarkonium provides an ideal 
laboratory for testing the perturbative and nonperturbative effects of QCD. An effective 
theory suitable for describing these systems is nonrelativistic QCD (NRQCD)[l|, which is 
derived from QCD by considering the underlying nonrelativistic properties. According to 
NRQCD factorization fl, decays of heavy quarkonium into light hadrons or photons can be 
organized in a hierarchy of long-distance matrix elements (LDMEs), which are classified in 
terms of f , the relative velocity of the heavy quarks in heavy quarkonium. 

Decays of ^5*0 heavy quarkonium into light hadrons (LH) and two photons are among the 
simplest processes. The short-distance coefficients for corresponding LDMEs at leading order 



in V have been computed previously to next-to-leading order (NLO) in 



4l2l|. Moreover, 



that coefficient for 77 decay has been calculated to next-to-next-to-leading order (NNLO) 



m a. 



13|. However, a.. 



to leading order in 



1 coefficients of LDMEs beyond leading order in v are known at best 



l6|,lli 



16(1 . It is well known that the calculation at leading order in 



suffers from large uncertainties due to strong renormalization scale dependence. Therefore, 
to give a more precise description for ^Sq heavy quarkonium decays beyond leading order in 
f , QCD corrections to these coefficients are apparently needed. 

In this paper, we will study QCD corrections to the coefficients of order v"^ LDME, 
namely, corrections at order ^ for ^5*0 quarkonium decays. Up to this order of corrections, 
there are two unknown LDMEs which should be fixed. Unfortunately, lattice calculation 
Cs l3, though based on first principles, suffers from large uncertainties. 



of these LDM: 
In Refs. 



19N2l[|. a new method was introduced to estimate LDMEs by combining potential 
models, lattice calculation, and experimental data. This method will also be used in this 
paper to determine the two unknown LDMEs. Then with our calculated a^f^ corrections 
we will be able to get updated estimates for the decay widths of ^Sq heavy quarkonium into 
light hadrons and two photons. 

The rest of this paper is organized as follows. We briefly introduce the theoretical pro- 
cedures for calculations of heavy quarkonium decays in Sec. [Tll In Sec. Illlt we describe 



kinematics and method of calculation for these processes. Results in perturbative QCD are 
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summarized in Sec. llVt while corresponding results in perturbative NRQCD are summarized 
in Sec. |Vl By using the matching condition, we give the updated short distance coefficients 
in SecJVT] to include our new agv"^ corrections. With these newly obtained results, we de- 
termine the two unknown LDMEs using potential models and make predictions for relevant 
decay widths in Sec. IVIII Finally, in Sec. IVIIIt we present a brief summary. 



II. DECAY OF HEAVY QUARKONIUM IN NRQCD 

The Lagrangian of NRQCD is derived from the QCD Lagrangian by integrating out the 
degrees of freedom of order mg, the mass of the heavy quark. Local 4-fermion operators are 
added to accommodate the inclusive annihilation decay of heavy quarkonium which happens 
at scale of order mg. The Lagrangian of NRQCD is 

-^NRQCD = -flight + -^heavy + 5C. (l) 

Here £heavy describes nonrelativistic heavy quarks and antiquarks and is given by 



X, (2) 



where ip is the Pauli spinor field that annihilates a heavy quark, x is the Pauli spinor field 
that creates a heavy antiquark and Dt and D are the time and space components of the 
gauge- covariant derivative D^. Terms corresponding to light quarks and gluons are given 
by flight and 

Aight = -^trG^^G^^ + (3) 

where G^y is the gluon field strength tensor, q is the Dirac spinor field for light quarks 
and the sum is over nj flavors of light quarks. Relativistic corrections to the basic effective 
lagrangian £heavy + -Ciight are included in 5C and its leading terms are those bilinear in the 
heavy quark or antiquark field, 

C2 



+ -\ (V^t(D . - • D)^ + x^D ■ (7E - (?E ■ D)x) 

+ UHiT> X gE-gEx iD) ■ aip + xH^D x gE - gE x iD) • ax) 

C4 



+ TT- (^^(^B ■ ct)^ - x^gB ■ a)x) , (4) 



where E' = G^' and B' = \e'^^&^ are the electric and magnetic components of the gluon 
field strength tensor G^'^ . 

Further corrections include the description of inclusive annihilation decay of heavy 
quarkonium and can be achieved by adding local 4-fermion interactions as 

-fermion V^^^^aM, (5) 

— TTLn 
n V 

where /ia is the NRQCD factorization scale, Onif^A) is the local 4-fermion operator, dn is the 
naive scaling dimension of the operator and /n(/iA) is the short-distance coefficient which 
can be calculated perturbatively. 

Thus the decay width of heavy quarkonium can be given by the following factorization 
formula 

r{H) = J2?^^{H\OM\H), (6) 

n V [—1 

where heavy quarkonium state in the Fock space can be written as [2j 

iHC'^'Lj)) = 0{1)\QQC'^'L^)) 

+ Oiv)\QQC'^'L±l%) 

+ Oiv')\QQr'^'L%) (7) 

+ 0{v')\QQC'^'L^j'%g) 

+ 

and the relative importance of the 4-fermion operators regarding v can be accessed through 
the velocity scaling rules outlined in Ref.j2|. We conform to this standard NRQCD power 
counting rules throughout this work, although alternative power counting rules exist 
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24| . A detailed discussion of the influence of different power counting rules can be found in 



Ref. 20|. For ^5*0 heavy quarkonium decays at order v"^, we need only consider the dominant 



^5*0 Fock state and two singlet operators with dimension 6 and 8: 



OCsi!^) = i^^xH. (8a) 



for light hadron decay, and 



^('^i'') = ^ [v^W(-tD)V + h.c. , (8b) 



Oem('5™) = V^^x|0)(0|xV, (9a) 



1 



V^^x|0)(0|x^(-fD)2^ + h.c. , (9b) 



for electromagnetic decay. For a generic color-singlet operator of the form On = 
ijj^ JC'^XX^ ^n'^ , applying the Vacuum- Saturation Approximation j2|, we get 

{H\On\H) = J2{H\^^'^nX\X){X\x^lCnm) 
X 

^ {H\^|J^ICUm{0\x^lCnm), (10) 

where the omitted terms are of relative order f ^ and are irrelevant of our calculations here. 
Therefore, in the following we use the notations 

{OCs;1^))h := {HCSI!^)\0CSI!^)\HCSI!^)) ^ {HCSI!^)\0M'SI;^)\HCSI!^)), (11a) 

{VCS^^^))h := {HCS^^^)\VCS^^^)\HCS^^^)) ^ (if(i5W)|PEM(^5W)l^('5W))- (lib) 

Then the decay width of ^5*0 heavy quarkonium at order f ^ is 

mCsH^) LH) = ^^{OCsS'))n + ^^(^(^4^^))^, (12a) 

mCsi:^) ^ 77) = ^^^^^{OCsI:))h + (12b) 
where the leading order LDME is related to the wave function at the origin as 

{OCS^^^))h = ^|7^^.(0)|^ (13) 
and a definition of the ratio of LDMEs is important in this work [20] 

mm) • 

where q is half the relative momentum of the heavy quark and anti-quark and it is also 
convenient to define 

{v'Ih = {q'ln/ml^. (15) 
To calculate the short distance coefficients F and G in Eq. ffT2l) . we use the matching method 

n 

|2|. Since the short distance coefficients are insensitive to the long distance dynamics, we 
can substitute the bound state with a pair of on shell quark and antiquark separated by a 
small relative momentum and exploit the equivalence of perturbative QCD and perturbative 
NRQCD to determine the short-distance coefficients 

AiQQ^QQ) = J]^^^(QQia(/^A)|gQ) ■ (16) 

pert QCD mf^ pert NRQCD 

n V 

The left side of this matching equation can be calculated pertubatively in QCD and the 
right side can be calculated perturbatively in NRQCD. Then, we can get the short distance 
coefficients fnifJ^A) whose imaginary part gives F and G in Eq. (|T2|) . 



III. KINEMATICS AND METHOD OF CALCULATION 



We work in the rest frame of the heavy quarkonium and assume the following notations 
for the momenta of heavy quark and antiquark 



Pq -- 
Pq 



2^ + ^' 



where 



P = (2Eq,0), 

q = (0,q), 



(17a) 
(17b) 

(18a) 
(18b) 



and En 



In our calculation, we adopt the covariant spin-projector method [2514271] to project out 



the spin-singlet amplitudes. The projector we use is 



27| 



2v/2(E, + .n,)^2+^ + "^«) 8El ^2 



(19) 



To expand the decay width in terms of q, we make the following rescaling for any momentum 

k ^ k'Ejmg, (20) 

which leads all momenta independent of q, that is, dk'- ■ k'j/dq = 0. Thus we can expand the 
amplitudes in q before loop integration and phase space integration and extract the S-wave 
contribution by making the replacement 



q^ 



p'p' 



(21) 



D-1' ' Ami'' 

where P'^ is the rescaled momentum of the heavy quarkonium which equals (2mQ, 0) in its 
rest frame. Contributions coming from potential regions in perturbative QCD and pertur- 
bative NRQCD cancel each other exactly so we neglect these terms to simplify calculations. 



IV. PERTURBATIVE QCD RESULTS 



We use FeynArts [28|, l29| to generate Feynman diagrams and amplitudes and use self- 



written Mathematica codes to perform the remained calculations. Ultraviolet and infrared 
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f 



f 




FIG. 1: Feynman diagrams for ^Sg^' — )• gg at Born level. 

divergences are regularized with dimensional regularization and D = 4 — 2e is assumed. 
Ultraviolet divergences are removed by renormalization. We define the renormalized heavy 
quark mass mq, heavy quark field tpq and gluon field in the on-mass-shell shceme(OS) 
and define the QCD coupling constant g in the MS scheme, that is, 



9' = Zfg, ml = ZZlmQ, = Z^^^q, AI=^Z0^A„ (22) 
where terms with superscript denote bare quantities and Zi = 1 + 6Zi with SZi given by 

(23a) 



SZ^l = -3C^f^/,[— + | + 21n(2)]. 
An euv <j 

5Z^^ = -CF^M— + — + 6\ni2)+A\, (23b) 
4vr euv em 



SZ. 



OS 



Air 



f,mnf)-2CA)i 



euv eiR 

2 



SZ. 



MS 



•f^^/.l^ + ln(=|) + 2M2)l, 

47r 2 euv l^f- 



(23c) 
(23d) 



where = r(l -|- e)[^^^]% (3o{nf) = yC^ — ^TpUf, Hr is the renormalization scale and 
n/ is the number of light quarks. 



LH 



At leading order in as, there are two diagrams as shown in Fig. [H The corresponding 
Born level decay width and its relativistic corrections are 



2! 2(2mQ) 



{2D) 



>.4vr)' 



16 
9mQ 



;i-2e)(l-e)(0(i4^1))Lo, (24a) 

(24b) 



ri(^4'U^^) = -^^T^CS^^^^gg). 
3 m^, 



(2D) = 'i~ r{2-2 \ two-body phase space for q = in D dimension. 

Baa, 



where $ 

Our results agree with those in Refs 



9|. At next to leading order in a^, there are 
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{ 



( 



FIG. 2: Representative Feynman diagrams for —J- gg at one-loop level. 



virtual corrections and real corrections. Fig. |2] corresponds to Feynman diagrams of virtual 
corrections, where only distinct forms of diagrams are shown. Contribution of these virtual 
corrections reads 



^^r.CsfU,,)/.{|-i + (-il + !^)i 

+ —(-44 + 197r2 + {Anf~ 66) ln(^))] 
36 fif. 

q\4 1 4 , 44 , , , 



(25) 



1 n^'n n 

+ —(-4(11 + 241n(2))n/ + 24(33 - 2n/) In(^) - 26 77r2 + 874)]}. 



324 



While other terms agree with those in Refs.j?, 9|, the result of relativistic correction here is 
new. Feynman diagrams for real corrections are drawn in Fig. |3] and Fig. HI which correspond 
to final states with three gluons and qqg respectively. Results for these two sets of real 
corrections are 



-riicK . \ CaOs , 1 11 181 23 2 

^{ S'o' ^ ggg) = /erB( S'o ^ gg){- + — + — - —tt 

71 6e 18 24 



(26a) 



Ti^S"^^ ^ qqg) = UfV^i'S','' ^ gg) 



e 54 
a 



fe 



2 16 q2 8 104 

TT rri + e)r(i ^ ^^'s'e ~ Y ^ ^^9^ ^ 'W^'' 



(26b) 

While other terms agree with those in Refs.j?, 9|, the results of relativistic corrections here 
are new. Adding Eqs. ( |25l) and ( |26l) . we obtain the NLO QCD corrections plus relativistic 
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f 



f 



FIG. 3: Representative Feynman diagrams for ^S^^''^ — t- ggg. 



FIG. 4: Representative Feynman diagrams for ^ qqg 



corrections for the light hadron decay width of ^Sq heavy quarkonium 



■^NLO {1 q[1] 

qcdI 



S'o ^ LH) 



216' 



-64n/ + 12(2n/ - 33) ln(^^ 



I 1 nnol I ^ r I ^nA^^^t '-^ 



- 937r^ + 1908] + ^^[^ + ^(24(ln(-^) + 21n(2))(33 - 2nj) (27) 
+ 164n/ + 23771^-4964)]}. 

Adding these terms together, we get the ^Sq decay width into hght hadrons in perturbative 
QCD 



TqcdI^^W ^ LH) = rB(^4'^ ^ gg) + fK^^^J ^ ^^) + T^^^cbC S^^^ ^ LH). (28) 



For QCD corrections to the electromagnetic decay, there is no real correction. Diagrams 
at Born level and one-loop level are the same as those in Fig. [1] and Fig |2] except that 
the final state gluons are substituted with photons and diagrams containing three-gluon or 
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four-gluon vertexes are excluded. We then get the results 

r|(^4'' ^ 77) = -|4-rB(^5W ^ 77), (29b) 

^ /I ^fil N ^ /I ^fil N . «s rTT^ — 20 ,16 196 — ISvr^,, 

rv('4 ^ ^ 77) = rB(^sW 77)/.^[^— + + ^7 )], (29c) 

where eg is the electric charge of the heavy quark. Results of relativistic corrections in Eq. 
f l29cp are new and the other results agree with those previously calculated as summarized 
in Ref. 0. Adding Eqs.f l29a| ). ( I29b[) . and f l29cp . we get the result for the 77 decay width of 
^5*0 heavy quarkonium in perturbative QCD 

rQCD(^4'^ ^ 77) = rB(^4" ^ 77) + rK^^w ^ 77) + rv(^4" 77). (so) 



V. PERTURBATIVE NRQCD RESULTS 

Order a^f^ corrections to the leading order LDME ((9°(^S'q^')) in perturbative NRQCD 
have been calculated in Ref. j2|, where a cutoff was introduced to regularize the ultraviolet 
divergences. We rewrite it in dimensional regularization, 



6iT /i^ euv eiR rrif- 



We define the renormalized operator O^i^sl^^) using the MS scheme 

C»°(i5W) = zgfSoR(i^W), (32) 



where 



^ Ft + ln47r - 7E 

Stt Hi euv ml 



Z^' = l- Ii;£r^(- )^(_ + in4vr - ^^)^. (33) 



Therefore 



{0YSS%L0 = [1 + ^{\n- + ln4vr - 7e)\](0(^4^'))lo. (34) 
Considering that 

{VCS^^)ho = q^(0(^sW))LO, (35) 
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the decay width into hght hadrons in perturbative NRQCD becomes 



rNRQCD("'^'S'o^^ — > LH) — 

{FCSS:') + ^[GCS^^') + '-^4n- + ln4vr - 7E)F(^5W]}i^i^4^- (36) 

The electromagnetic decay rate can be obtained by replacing F{^Sq^) and G'(^5'o^') with 
F^^{^Sq^) and G^^{^Sq^) respectively. 



VI. MATCHING 

Finally we obtain the short distance coefficients by equating results from perturbative 
QCD in Eqs.([2H]) and ^ with that from perturbative NRQCD in Eq. 

m 47ra2 -64nf + 12(2nf-33)ln(^)- 937r2 + 1908 
^(i^W) ^^^^ ^ ^ i ], (37a) 

GCS^^^) =^{-l + ^^[481n(2)(25 - 2n,) + 164n, - 4964 

2 2 

+ 24(33 - 2nf) In(^) + 192 In(^) + 237n^]}, (37b) 

F, ,CS^^^) =2vra^eJ(l + ^^^), (37c) 

G, ,(1SW) =27ra%{-^ + ^l[481n(4-) - 961n(2) - ISvr^ + 196]}, (37d) 

where QCD corrections for Gf^Sn^') and GyyCsj^^) are new while the other results agree 

n n p n 

with those in Refs. [2|, |6|, 17|, 19[ . With these short distance coefficients, we can update the 
decay widths of ^5*0 heavy quarkonium into light hadrons and two photons. 



VII. PHENOMENOLOGY 

The above obtained result can be used in ^5*0 charmonium and bottomonium decays. In 
the following we will focus on the rjc decay width into light hadrons (approximately the total 
width) and decay width into two photons. In these decays there are two unknown LDMEs. 



In principle, one can fix these LDM 



calculation from lattice QCD 



17 



]s either through direct fit with experimental data 30| or 



18| . The order f ^ LDME is ultraviolet divergent and needs 



to be regularized [20|. For lattice calculations this is performed by imposing a hard cut-off 

11 



regulator. However, due to slow convergence of this regularization, the results available 
from lattice calculations of order f ^ LDME suffer from large uncertainties 17|, [iSj . On the 
other hand, we find that direct fit of the two LDMEs using experimental measurements of 
77 width and total width of rj^ can not give reliable values due to the approximate linear 
dependence of the two theoretical predictions for these two decays. 

Therefore we determine the two LDMEs using the potential model method recently in- 



troduced in Refs. 



19l-|2l|. A widely accepted potential model, the Cornell potential jsil 

V{r) = -- + ar, (38) 



is chosen in this work. Since the spin dependent effect is not included in this potential, the 
LDMEs calculated this way are accurate up to corrections of relative order v"^. However, as 
argued in Ref. [20^, this error is in fact much less than the order (about 30%), thus we 
attach an uncertainty of 30% to the central value of the order f ^ LDMEs to account for the 
error due to this static potential approximation. 



In solving the Schrodinger equation 32|, there are three unknown parameters, a 
0.1682 ± 0.0053 GeV^ is taken from the average of lattice calculations j20| and the mass 

2o|. Here we take m(V' (25))- 



19 



parameter is expressed in terms of the 15-25' mass splitting 

m{J/ilj) = 589.188 ± 0.028 MeV 33|]. The last remaining parameter is fixed by equating 
theoretical predictions to experimentally measured results. When we use the decay width 
formula, we resum a class of relativistic corrections at leading order in for 77 decay as in 
Refs. 20|, |2l|]. For the experimental input, we make use of this approximation r^^{ric{nS)) = 
r*°*'^'(?7c(?T.5)). For rjc, we use V^'^lrj^) (or r^^(?7c)) as input to obtain one set LDMEs which 
are then utilized to obtain T^^^rjc) (or V^'^^rjc)). For ric{2S), we use the total width as input 
and make predictions for the 77 decay width. We take nic to be 1.4±0.2 GeV |20|, a = 1/137, 
Aqcd = 0.39 GeV and vary the renormalization scale fir and NRQCD factorization scale 
/iA separately from 1 GeV to 3 GeV with 2 GeV as the central value. The LDMEs are 
expressed in terms of the wave function at the origin |7^(0)p and {v"^) as defined in Eqs. 
( !T3|) and (fT5|) . In each determination of these LDMEs and the corresponding decay width, 
we evaluate the variations caused by the uncertainties of the parameters and summarize 
them in the tables. Because the potential does not take into account of the spin effects, we 
attach each LDME an uncertainty (i?^) of 30% of the central value. In each case, various 
uncertainties are added in quadrature to give the total uncertainty. 
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For rjc, with the 77 width T'^'^{ric) = 7.2 ± 0.7 ± 2.0 KeV [33] as input, the determined 
LDMEs are 

|7^;^(0)r = O.SSlt'oill GeV^ (39a) 

{vX = 0.2281°;}^^, (39b) 

where the superscript 77 indicates that we use the 77 decay width as input. Various 
uncertainties are summarized in Table [H The most significant uncertainty comes from the 
experimental data. At order a^f^, another dependence on NRQCD factorization scale fi\ is 
introduced. However, as we can see from Table [H variations of LDMEs are small when we 
change from 1 GeV to 3 GeV. In Fig. O we present fir dependence of three sets of LDMEs 
in terms of |7?.^;,(0)p and {v'^)^^^, where the other parameters are fixed to their central values. 
Of these three sets of lines, LO represents calculation without any QCD corrections, NLO* 
corresponds to that including QCD corrections but only for terms at leading order in v, 
and NLO means our new result with order a^f ^ correction taken into account. The LDMEs 



corresponding to NLO* have been computed earlier in Ref. 20| and can be compared here 
with the values including the new order a^f ^ corrections. We can see from this figure that 
these two lines are close to each other, which refiects the fact that the effect of the new order 
^ correction is not large. Utilizing this set of LDMEs as input, we get the total decay 
width for rjc 

ptotai^^j = 31 Atlll MeV. (40) 



This value is in consistency with experimental measurement 28.6 ± 2.2 MeV [331], although 
there are large uncertainties. The details of the uncertainties are summarized in Table [B 
Here, the uncertainty induced by the fir dependence predominates over that from the ex- 
perimental input of the 77 width of rjc. 



If, on the other hand, we use the total width of r/^, T^"^''\r]^) = 28.6 ± 2.2 MeV [331 as 
input, then we get another set of LDMEs 

|7^Jf(0)p = 0.814t°ii GeV^ (41a) 
{vX^ = 0.23Atroll (41b) 

where the superscript LH indicates that we use the total width of rjc as input. Variations 
with respect to the parameters are summarized in Table [Tll The experimental uncertainty in 
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(a.) fir dependence of |7^Jl|J(0)p 
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0.24 - 



0.22 



0.20 




(b)/ir dependence of 

FIG. 5: fir dependence of LDMEs for rfc using the observed 77 width as input. LO represents 
values without QCD corrections, NLO* includes QCD corrections only for terms at leading order 
in V, and NLO takes into account our new QCD corrections to order v'^ terms. 

this case is small and the main uncertainty comes from the relatively strong Hr dependence 
of the theoretical prediction. The fir dependence of the two LDMEs is shown in Fig. [61 As 
in previous case, we display another two sets of LDMEs, where only QCD corrections at 
leading order in v are taken into account or no QCD correction is considered. The two sets 
with QCD corrections show great improvement of fir dependence with respect to the one 
without QCD corrections, and they are almost parallel to each other. The only difference 
between the two sets of values is that the a^f^ correction enhances 17^7)^(0) p by about 30%. 
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TABLE I: LDMEs obtained from potential model for r]c and the predicted total width of r/c, the 
superscript 77 indicates that the observed width of rjc — )• 77 is used as input. The second row 
gives central values while the followed rows give variations with respect to related parameters. 
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This enhanced LDME coincides with previously obtained value using T^'^{ri^ as input. The 
value of {v'^)r^^ is relatively stable. With this set of LDMEs, we obtain the 77 decay width 



r^^(r^,) = e.eilLl KeV, (42) 



33|. 



which is also consistent with the experimental measurement 7.2 ± 0.7 ± 2.0 KeV 

Since now we have two sets of values for the two LDMEs for rj^ in Eqs. (l39|) and (l^TI) . 
we can combine these values to get a better estimation. The uncertaintes in Table [T] and 



Table [TTl are correlated, and we use the method in Ref . 201] to treat these correlations. First 
we construct a two-by-two covariance matrix for {Oi^SQ^))^^^ {0{}S^q^))]:^. It describes 
correlations between the variations in the two tables and is defined as Cjk = Xlt ^ji^ki with 
Aji = |(0^ — Ojj). The indexes j,k refer to these two leading order LDMEs and i runs 
through every item in Table [T] and Table [TTl Oj] and Oj^ correspond to the plus and minus 
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|«LH(0)|2(GeV') 




(a.) fir dependence of |7?.)^^(0)| 




{h)fir dependence of (v^)^^ 

FIG. 6: fir dependence of LDMEs for r/c using the observed total width as input. LO represents 
values without QCD corrections, NLO* includes QCD corrections only for terms at leading order 
in V, and NLO takes into account our new QCD corrections to order v"^ terms. 



variations of the LDMEs. For the i-th item in Table [T] or Table [TTl we define the xl &s 



0i){C-^)Jk{0k^-0i 



(43) 



and minimize it to get the average value Oj for {p(}S^Q^)n^. Once we obtain the values of 
{0(}S^Q^)rf^, we use the potential model to get the values of {V{}S^Q^)r,^. We perform this 
calculation for each of the items in Table [T] and Table |Tl] , treat the renormalization scale 
/i^ and NRQCD factorization scale /ia simply as the same quantities in the two tables and 
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TABLE II: LDMEs obtained from potential model for r]c and the predicted 77 width of r]^, where 
the superscript LH indicates that the observed total width of r/c is used as input. The second row 
gives the central values while the followed rows give variations with respect to related parameters. 
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4- A /T 




01 


1 Q 


-Act 


-0.038 


-0.018 


-0.16 


+A/i,. 


0.201 


-0.017 


2.47 


— Afir 


-0.189 


0.016 


-2.73 


+A/iA 


-0.020 


0.002 


0.13 


-A/iA 


0.036 


-0.003 


-0.21 


+Ar,, 


0.052 


-0.004 


0.46 


-Ar,, 


-0.053 


0.005 


-0.47 



express the results in terms of |7?.,,^(0)| and (f^),,^. The results are 

i^..(o)r 



0.834t°:?«i GeV^, 



0.232 



h0.121 
-0.098' 



(44a) 
(44b) 



with the details of the uncertainties given in Table IIIII We note that the uncertainties here 
for |7^r,^(0)p are smaller than those in Eqs. (15^ and ( HTj) . 



For ric{2S), we use the observed total width 14 ± 7 Mev|33| as input and get the LDMEs 
for T],{2S) 



17? 

\'Sci2S) 



(0)1 



.2\LH 



Q ociCi+O.lSO 



(45a) 
(45b) 



Some potential model calculations of the squared wave function at the origin in Ref . |34j| give 
0.418 for the logarithmic potential[35|] , 0.529 for the QCD-motivated B-T model sj, and 
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TABLE III: The averages of the LDMEs for r]c- The second row gives the central values and 
subsequent rows give variations with respect to various uncertainties. 



Case 


\n„ (0)1 2 (GeV^) 

1 '/c V /l J 


\ / '/c 


central 


0.834 


0.232 


' \ / ijc 


0.159 


0.070 


\ 1 '{C 


-0.115 


-0.070 


+Amc 


0.170 


-0.065 


—Arric 


-0.113 


0.097 


+Ao- 


0.036 


0.019 


-Ao- 


-0.033 


-0.018 


+A/ij. 


0.124 


-0.010 


— AUr 
' — 


-0.069 


0.006 


+A/iA 


-0.025 


0.002 


-A/iA 


0.044 


-0.004 




0.069 


-0.006 


-Ar;^: 


-0.072 


0.006 




0.037 


-0.003 


-Ar*"*'^' 


-0.037 


0.003 



0.559 for the power-law potential js^]. Our result is consistent with their values. Table HVl 
gives details with various uncertainties. In the ric{2S) case, both the /i,. dependence of 
theoretical result and the experimental input of the total width have large uncertainties and 
therefore the LDMEs are subject to relatively large uncertainties. In Fig. [TJ we present fir 
dependence of this set of LDMEs. The shape of the lines is similar to Fig. [6] except the value 
for |7?.^^(25)(0)p here is smaller by about a factor of 2. Exploiting this set of LDMEs, we 
can make predictions for the 77 decay width of r]c{2S), 

r^^(77,(25)) = SMtlfo KeV. (46) 

This prediction is consistent with the experimental observation that the branching fraction of 



77 decay is smaller than 5x10 ^ [33|. Another experiment measured T^^{ric{2S))B{ric{2S) — )■ 



KKn) = (0.18 ± 0.05 ± 0.02)T^^{r]c{lS))B{r]c{lS) KKn 



and assumed that the 
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TABLE IV: LDMEs obtained from potential model for rjc{2S) and the predicted 77 width of 
r]c{2S), where the superscript LH indicates that the total width of r/c(25) is used as input. The 
second row gives the central values while the followed rows give variations with respect to related 
parameters. 



Case 


|7^LH(o)|2 (GeV3) 




r(r/,(25) ^ 77) (KeV) 


central 


0.423 


0.255 


3.34 




0.121 


0.076 


0.58 




-0.077 


-0.076 


-0.33 


+Amc 


0.077 
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-0.08 


—Anic 


-0.035 


0.099 


0.47 
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-0.019 
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+AfJ,r 
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-0.016 


1.22 


-A fir 


-0.091 


0.014 


-1.33 


+A/iA 


-0.012 


0.002 


0.07 


-A/iA 


0.022 


-0.003 


-0.11 




0.167 


-0.026 


1.47 


-'^rr,,(25) 


-0.192 


0.029 


-1.58 



branching fractions of rjc and ric{2S) decays into KsKn were equal, and made use of T{ric — >■ 
77) = 7.4 ± 0.4 ± 2.3 KeV, and then derived T{ri^{2S)) = 1.3 ± 0.6 KevU, y. Our result 
is not in contradiction with their measurement within errors. 



VIII. SUMMARY 

Within the framework of NRQCD, we calculate order a^f^ corrections to decays of ^Sq 
heavy quarkonium into light hadrons and two photons. In both processes, infrared diver- 
gences are found to be canceled through the matching of perturbative QCD and pertur- 
bative NRQCD results. There are two unknown NRQCD LDMEs, which are determined 



using potential model method 19l-l2l| either with the observed total width or two photon 
width as input. When using r^^(r/c) as input, we get |7^;57(0)|2 = 0.88ll°:i^®^ GeV^ and 
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|«;;«J0)|2(GeV') 

0.9 
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^,(GeV) 

(a)^r dependence of |7?.5^;^2S) (O)l^ 
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(b)^r dependence of (■y^)Jj^2S) 

FIG. 7: dependence of LDMEs for rjd^S) using the total width as input. LO represents values 
without QCD corrections, NLO* includes QCD corrections only for terms at leading order in v, 
and NLO takes into account our new QCD corrections to order v'^ terms. 

{v^)^'] = 0.228l|J:i^[;, from which we predict r*°*^'(r/,) = 31.4+^9-3 MeV. Alternatively, when 
using r*°*^i(r/c) as input, we get |7^^f (0)|2 = O.SUtoHl GeV^ and {v^)l;^ = 0.234+°;J2i, and 
we predict the 77 width of r]c to be T'^'^('r]c) = 6.6II2 33 KeV. All these predictions agree well 
with experimental data. We then combine these two kinds of determination of LDMEs and 
get the average values |7e^,(0)|2 = OMAtoHl GeV^ and (v^)^^ = 0.232l:oml- For r]c{2S), 
we use the observed total width as input and find |7^|^^^2S)(0)P = 0-423+o;23o GeV^ and 
(^^)S2S) = 0.255;o:}o9. With this set of LDMEs, we predict the 77 width of ?7c(25) to be 
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3.34l2'io KeV, which is not in contradiction with data within uncertainties. Consequently, 



-2.10 

,2 



the order agV corrections (especially the one to the decay into light hadrons) are found 
to have significant effects on improving the consistency between theoretical predictions and 
experimental measurements. 

Note added. When we finished the calculations and are preparing this paper, a related 
work appears |39| that also gives a^w^ corrections to the 77 decay width. We find our results 
for this channel agree with theirs, while we have also calculated the light hadron decay 
width. 
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